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In a general framework that has been labeled the "gauging of equations method", we study the 
diagrams that contribute to Compton scattering off a relativistic composite system. These con- 
tributions can be derived for iV-particle bound states described by the covariant Bethe-Salpeter 
, equation with a method equivalent to minimal substitution in the one-particle case and yield the 

correct contributions (including subtraction terms) in the order 0(e 2 ). We give the Ward-Takahashi 
identities for the general two-photon vertex as well as the corresponding constraints for the two- 
photon irreducible interaction kernel and the Bethe-Salpeter amplitude describing the bound state. 
From this we can show that gauge invariance holds for the full two-photon vertex. We furthermore 
ON I study in detail the low-energy limit of the Compton scattering tensor in this approach (including a 

discussion of the pole terms) and can prove that the full amplitude yields the correct Born-Thomson 
limit as we shall explicitly show for the spin-0 case. 

The calculations are completed by the investigation of certain approximations that can be formu- 
lated for arbitrary TV-particle bound states. We neglect for instance contributions from n-photon 
irreducible interaction kernels and show that in this case gauge invariance is only realized if either 
, the interaction kernel in the Bethe-Salpeter equation is independent of the total momentum and 

■ additionally is of local type, or if the photon energies vanish; furthermore, we find the correct low- 
energy limit in this approximation. To clarify our approach, we also give the results in the order 
0(e); as examples, we will quote some resulting lowest order expressions for a qq system explicitly. 

o 

■ I. INTRODUCTION 

The first order result for a Green's function in an external field obtained by a minimal substitution prescription and 
the analoguous first order result for Bethe-Salpeter amplitudes have been shown in refs. jl]-|| by Kvinikhidze and 
Blankleider; the authors refer to this procedure as "gauging a hadronic system" . However, up to now nothing is known 
. [ about the important question how this general scheme can be applied to Compton scattering off a hadronic system 
which is a second order process. In fact, this is a difficult task since electromagnetic gauge invariance is intimately 
related to the dynamics of the strong interaction inside extended composite hadrons (see [|J and references therein). 
This important question has been addressed by Ito and Gross in ref. [pd[| where the authors have derived various 
diagrams that have to be included to guarantee a gauge invariant description of the Compton scattering process. 
These considerations — based on a diagram analysis in a general four-dimensional Bethe-Salpeter framework - were 
explicitly applied only to the deuteron system (i.e. with only one charged constituent) as an example that keeps 
the equations as simple as possible. In contrast to this approach, the procedure of gauging the relevant bound state 
equations as applied here has the advantage that it allows for a complete overview of the possible contributions 
in a given order 0(e n ) for any TV-particle bound state; it thus provides a more systematic view on the numerous 
processes contributing to Compton scattering off a bound state. In this way, it is therefore possible to study certain 
approximations (like e.g. neglecting n-photon irreducible interaction kernels) on the basis of the knowledge of the 
full expressions; this seems to be a more firm way of introducing and handling these approximations than in the case 
where one does not know what terms are precisely suppressed under certain assumptions. 

In this article, we aim — on a general level — at a complete classification of the diagrams that contribute to the 
Compton scattering process A'f — > Aj for bound states A irrespective of their composition. Although our approach 
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can be considered as a quite technical study of this process, it seems worthwhile to present it in this publication since 
it provides a more transparent way of discussing this issue compared to the technically involved explicit calculations 
in ref. pH ; furthermore, it allows for a simple check of the correct low-energy limits. We thus demonstrate that the 
method of Blankleider and Kvinikhidze (see Jl]-|| ) turns out to be a powerful tool to derive all contributions not only 
in first order of the electromagnetic interaction but also in the order 0(e 2 ). Herein, the relevant degrees of freedom 
are the constituent particles of the composite system; our results for the full two-photon vertices turn out to be 
independent of the interaction kernel used in the underlying Bethe-Salpeter equation describing the bound state. In 
order to make clear our notation and the scheme that we apply to the second order amplitudes, we will always give the 
results in first order (i.e. 0(e)) for comparison. To simplify the notation, we will mostly suppress the dependencies 
of relative momenta between the constituents; however, we will give some explicit examples for qq systems including 
the full dependency of the internal relative momentum. 

After some introductory remarks on fundamental equations, definitions and notations in section [il], we will derive 
the second order expansions of Green's functions, two-photon vertices and Bethe-Salpeter amplitudes in section III. 
Ward-Takahashi identities and differential expressions with respect to the bound state's total momentum are given 
in section ^ for vertices and n-photon irreducible interaction kernels. The issue of gauge invariance is addr esse d in 
section ^ before we study the low-energy limit of the Compton scattering tensor in section VI . In section VII , we 
summarize our results. 



II. GREEN'S FUNCTIONS, T-MATRIX AND BETHE-SALPETER EQUATIONS 



The fundamental Green's function equation for a given number of n fermions and n anti-fermions — depicted in fig. 
(0) in a graphical form — is defined as follows: 



G = Gq — iGqKG 
= Go — iGKGo 



(1) 



G 







; k Ec 


G 







FIG. 1. The equation for the Green's function G according to eq. (|l|). The vertical dots (here and in the following figures) 
denote the possible occurence of additional fermion lines depending on the number of constituents. 



Throughout this article, all n + n — 1 four-dimensional integrations over internal relative momenta are implicitly 
understood; furthermore, we omit all indices for internal quantum numbers such as flavour, colour and spin. To clarify 
this point, we write down the fundamental Green's equation, i.e. the first line in cq. ([!]), for G (P; {kj}; {kj}) = 
G {{p'i}\ {pi}) with i = 1, 2, N and j = 1, 2, N - 1 explicitly as 



[G (P; {*<}; {A*})] {aM} = [Go (P; {kfr {%})] { ^ } 



(2) 



(2tt 



,4(JV-1) 



(27T) 



4(JV-1) 



[Go (P; {fc;}; {g}})] {Qi}>{74} [K (P; {<£}; {q })] ^ [G (P; {q,}; {M)W, {ft 



where we have separated the total momentum P = pi + . . . + Pn = p'\ + • • • + p'n {N = n + n) and we integrate over 
all relative variables {kj} = k\, k%, . . . , fcjv-i; see also appendix [A|. Here, the index set {a>i} = ai, . . . , ajy represents 
the multi-indices for each constituent particle in flavour, colour and spin space. Since Go describes free propagation 
without any interaction between the constituents, it is defined as 



G (P; {kfr {kj}) := G (P; {%}) • (2n) 4 6* (k{ 



(27T)- 



S 4 (k' N ^ 



in this explicit formula. Here and in what follows in this section, all quantities are functions of the total momentum 
P; furthermore, the dependence on the relative momenta kj of the constituents (aside from some examples) and all 
indices will be suppressed in this paper. 

The kernel K contains all irreducible interaction terms between fermions and anti-fermions, se e fig. (|^). In this 
publication, we will not make any assumption about this kernel; however, we will show in section IV B| that it has to 
obey certain - general — Ward-Takahashi identities depending on the approximations applied in order to preserve 
gauge invariance for electromagnetic processes. 
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FIG. 2. Irreducible diagrams that contribute to the interaction kernel GoKGo with N — 4 constituents as an example. 



Let us stress that this kernel K includes not only iV-particle irreducible diagrams but also diagrams in which some of 
the constituents are unaffected by the interactions between the other constituents; an example is given in the second 
diagram in fig. (|J). In this sense, the kernel K in the Green's function equation is the sum of all possible i-particle 
irreducible kernels {i = 2, . . . , N) (see e.g. fit] for a precise formulation of K = ^T ifW in the case of qqq baryons). 
The iteration of an infinite sum of irreducible diagrams in eq. (|l|) can also be written as G = Go — iT by introducing 
the T matrix as T := GqKG which then satisfies the following equations: 



T = GqKGq - iG KT 
= GqKGq - iTKGo 



(3) 




T 11= 1 K 



= 




FIG. 3. The equation for the T matrix in eq. (E 



From the equation for the Green's function, we can deduce the Bethe-Salpeter equation for the Bethe-Salpeter 
amplitude x describing on-shell bound states and its adjoint x (see p^Jl^ ]): 



X = -iG Kx and x = ~ixKG a 



(4) 



Here, x is a short-hand notation for the Bethe-Salpeter amplitude Xa 1 ,...,a n+R {P', {kj}) that in coordinate space is 
defined explicitly as 



Xai,...,a„+B (>£lj • • • ! x n+n ):=(0 \ Ti> ai (x 1 )...i> an (x n )^ otn+1 (x n+1 )...^i 3n+fi (x n+ f i )\ P) 



(5) 



e lPX I ^ e~ ik ^ . . . e- ik -^x a u..., an+ AP; {fc,-}) 



(27T) 



4(AT-1) 



for a bound state of n fcrmions and fi anti-fcrmions (N — n + n). In this definition, T denotes the time-ordered 
product; see also appendix |X| for the definition of the coordinates X, r*i, . . . , rjv-i- 

With the formal solution of eq. ([!]), i.e. G — (Gq 1 +iK the Bethe-Salpeter equations in eq. (Q) can be formulated 
as follows: 



G^x = and xG" 







(6) 



Note that these homogeneous Bethe-Salpeter equations are only valid (and thus the Bethe-Salpeter amplitudes are 
only defined) for on-shell bound states of momentum P and mass M with the constraint P 2 = M 2 . 
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FIG. 4. The Bethe-Salpeter equation for a bound states of iV constituents in a graphical notation, see eq. m 
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III. SECOND ORDER EXPANSION IN THE ELECTROMAGNETIC FIELD 



In the following we extend the fundamental Green's equation and the Bethe-Salpeter equation by expanding operators 
= Go, G, K, x in the presence of an external electromagnetic field as 

A = - ieO 7 - e 2 0' n + 0(e 3 ) . (7) 

In this way, we find the expressions for the second order Green's function G 77 and the second order Bethe-Salpeter 
amplitude x 77 - 

A. Green's Functions 

The Green's function in eq. can be transformed into a Green's function in an external field via minimal substitution: 

G A = Gq — iG A K A G A (8) 

where eq. (^) for the operators = Gq,G,K is applied. In this way, each fermion line in the operators will be 
attached by a number of photons in the order 0(e n ) similar to the well-known construction of the fermion propagator 
S A (p,p') in an electromagnetic field, see fig. (||). 




FIG. 5. The fermion propagator in an external field. 



1. First Order Expansion 

With the expansion in eq. (^|) and by using eq. (|lj), we can write down an equation for the Green's function in first 
order of the electromagnetic coupling: 

G 7 = G 7 - iGlKG - iG K"<G - iG KG~< 
= (1 + iGoKy 1 (G 7 (1 - iKG) -iG K~<G) 

=GG" 1 =g- x g 
— G (G " 1 G 7 Gg 1 - iK 1 ) G . 

Defining the amputated full one-photon vertex T 7 := G _1 G 7 G _1 and the amputated free photon vertex Tq := 
G ~ 1 GqGq , we can write 

g 7 = Gr 7 G with r 7 = r 7 - ik* 1 (9) 

for the first order Green's function in an external field; this equation is depicted diagrammatically in fig. (^|). 




FIG. 6. The one-photon vertex GqF^Go according to eq. fflh in a graphical notation. 
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This full one-photon vertex can not solely be studied in models where the Bethe-Salpeter equation is used in its 
instantaneous approximation; it is also relevant in models including retardation effects (see S). The concept of an 
"interaction current" K 1 has been introduced in ref . jL7) ; its effects have been studied quantitatively in f| where the 
electromagnetic pion form factor is calculated with the full vertex T 1 in a four-dimensional separable ansatz for K 
based on 't Hooft's instanton-induced interaction (see also flic|). 

To make clear the physical meaning of K 1 , we give some diagrams contributing to this kernel explicitly in fig. (Q) as 
an example. Note that — in analogy to the kernel without electromagnetic field in fig. (|^) — all i-particle irreducible 
diagrams (i = 2,.,.,N) contribute to K 1 , i.e. also processes where some of the constituents can be regarded as 
spectators. 



K 1 — — 



+ 



+ 



+ 



FIG. 7. The one-photon irreducible interaction kernel GoK^Go with N = 4 charged constituents as an example. If the 
interaction kernel K includes charged exchange particles then additional diagrams contribute to K' 1 in which the photon also 
couples to these charged exchange particles. 



For bound states composed of more than two constituents {e.g. a nucleon described as a qqq system), the method 
of gauging the Green's function equation provides a solution to the notorious problem of overcounting diagrams in 
few-body systems. A detailed discussion of the role of subtraction terms emerging in this framework can be found in 
refs. HH where the problem of ambiguous cuts in the related diagrams (as well as its cure) is intensively studied. 
We remark moreover that our notion of a "lowest order vertex" with respect to Tq (and as well to T^ 7 introduced 
in the next paragraph) is somewhat sloppy since — apart from coupling constants that might be absorbed in the 
interaction kernels K, K 1 and K 11 — we have no ordering or power counting scheme that is reflected in this term. 
In principal, all parts of T 7 (and r 77 ) can contribute in the same magnitude to the full n-photon irreducible vertices. 
However, these lowest order vertices are minimal in the sense that they should be included in any framework that 
considers bound states like the pion or the nucleon as composite systems since they describe the simplest photon 
couplings to the constituents. Therefore, we will use this phrase now and then in the following discussion although it 
is not strictly exact. 



2. Second Order Expansion 

For the description of Compton scattering off a bound state described by the Bethe-Salpeter equation, we need the 
second order expansion of the Green's function. Let us therefore again start with eq. (|J) and evaluate it in the order 
0(e 2 ): 

G 77 = Gl 1 - iGj'KG - iGoK^G - iG^KG 11 - iG^K^G - iG^KG 1 - iGoK^C 
= (1 + iGoKy 1 (G 77 (1 - iKG) -iGoK^'G - iG^K^G - iG 7 KG~< - iGaK^G 1 ) 

=GG- 1 =G f 7 1 G 

= G{G^ 1 Gl 1 G^ 1 - iK~<"t -iC^ClE 1 - iG^G^KG^G" 1 - i/f 7 G 7 G~ 1 )G . 

The last three terms can be simplified with the expression for the full one-photon vertex in eq. (^): 
^77 = -iG^GoT^GoK^ - iG^GoT^ G KG T 7 GG" 1 - iK 1 GY 1 GG^ 1 

i(G-Go) 

= -ir^GoK-t + r 7 Gr 7 - r 7 G (r 7 - iK r ) - iK~<gt~< 
= r 7 Gr 7, -r 7 G r 7 . 

We define — in analogy to the first order expansion — the amputated full two-photon vertex r 77 := G _1 G 77 G _1 
and the amputated free two-photon vertex T 77 := Gq 1 Gq 1 Gq 1 ; the second order Green's function then reads 
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G 77 = gt 11 g with r 77 = r 77 - iK^ + r 7 Gr 7 - r 7 G r 7 



(10) 



and is shown in terms of diagrams in fig. (|8|) . This is the key result of this section since it provides — on the basis 
of rather simple considerations — the answer to the question which processes contribute to Compton scattering off a 
bound state. In this sense, eq. (|l^) can be regarded as a generalization of the results of ref. Jll]]; furthermore and in 
addition to it, it gives a more systematic overview and ordering of the various terms and by this means thus reveals 
the existence of a subtraction term besides the pole contributions in r 7 GT 7 . Note that this subtraction term T^GoT^ 
naturally emerges in this framework; it prevents overcounting of the lowest-order parts in r 7 GT 7 and turns out to 
be crucial for a proper description of Compton scattering off bound states. The two-photon irreducible interaction 
kernel K 11 in eq. ( |l0| ) can be understood analogously to fig. (jjj). 




FIG. 8. The two-photon vertex GqT~ i ''Gq according to eq. (jlOl) in a graphical notation. 



As an approximation of this full two-photon vertex, we will introduce G 77 := Gr 77 G where all n-photon irreducible 
interaction kernels, i.e. K 1 and K 11 , are neglected. Then, by virtue of G = Go — iT, the subtraction term can be 
evaluated explicitly — see fig. (||) — finally yielding 

g 77 = Gf 77 g with f 77 = r 77 - ir 7 rr 7 . (11) 




FIG. 9. The approximated two-photon vertex Gof 77 Go without inclusion of n-photon irreducible interaction kernels ac- 
cording to eq. (fill) in a graphical notation. 



3. Explicit Formulae for Compton Scattering 

These results in first and second order of the coupling e were formulated for a hadron in an amorphous electromagnetic 
field without stating precisely the photon momenta and coupling indices in Dirac space. If we consider Compton 
scattering kinematics (i.e. P + qi = P' + qz), we now have to indicate not only the correct momentum transfers but 
also to specify coupling indices /i, v, . . . of the photon lines and make sure that in each order all possible diagrams 
including permutations are taken into account. For the second order fermion propagator, this yields two diagrams 
(their sum being explicitly crossing symmetric) as shown in fig. (^0[). 




FIG. 10. The fermion propagator in an external field in second order with Compton scattering kinematics P' = P + qi — q-z 
for two photons 71 = 71(91, (J,) and 72 = 72(— q2, v). 
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For the first order Green's function in eq. (^|), all terms can be identified uniquely and therefore no correction is 
needed: 

G M (P', P) = G(P')T'"(P',P)G(P) (12) 
with V (P' ,P)=T% (P' ,P)~iK t " (P 1 , P) . 

Here we have introduced the momentum transfer P' = P + q for the photon 7 = 7(9, fj,) involved in this process. 
Considering the second order Green's function in eq. (fLCj), we find that additional terms appear due to crossing 



symmetry: 

G^(P', P) = G(P')T^(P', P)G(P) (13) 
with T^{P', P) = T% V {P', P) - iK" v (P', P) 

+ r»(P',P'- qi )G(P -q 2 )T v (P-q 2 ,P)-T%(P',P' - qi )G (P - q 2 )T%{P - q 2 , P) 
+ T V {P', P' + q 2 )G(P' + q 2 )T^(P + Ql ,P) - r^(P', P' + q 2 )G (P' + q 2 K(P + q U P) . 

For Compton scattering with the incoming photon 71 = 7i(<7i,m) and the outgoing photon 72 = r y 2 {—q 2 ,v), energy- 
momentum conservation requires P+qi = P'+q 2 where P (P') denotes the incoming (outgoing) total four-momentum 
of the bound state. After multiplying out all contributions from r 7 GT 7 , all terms in eq. ( p"3| ) can be identified with 
diagrams shown in |jllf where Compton scattering off a deuteron system (i.e. with only one charged constituent) has 
been studied. 

It will be helpful for some following remarks to split up the Compton vertex in eq. ( |l3| ) into parts that include the 
pole terms induced by the full Green's function G and the rest that is non-sing ular (see e.g. [||j24|): 



r<-(p', p) = r p ; LE (p', p) + r^ oPOLE (p', p) (w) 
with t^ ole (p',p) = t»(p',p' ~ qi )G(p -q 2 )v»{p-q 2 ,p) 

+ r"(P', P' + q 2 )G{P' + q 2 )T^{P + Ql ,P) . 

Let us note that from the equations above no constraint arises neither for the one-photon irreducible interaction 
kernel K fM (P',P) nor for the two-photon irreducible interaction kernel K^ V {P\ P). However, they have to obey 
certa in Ward-Takahashi identities in order to satisfy the constraints from gauge invariance as we will see in section 



1VB. Concerning the interpretation of the n-photon irreducible interaction kernels, it is worthwile to note that one 



can generally describe them in a different way than demonstrated above. Let us therefore first define 

G( e ) := G - it GqKG (15) 

as a modified Green's function that reduces to G( ) = Go describing free propagation for e = 0. With this definition, 
we introduce 



No 



= G„ 1 +ieK (16) 
= T%-ie K» 

= r(T - ie + r? e) G (6) rfo + r^G (e) i* e) - r^G r^ - rscrg 



and then finally 



y(e) 



= r (e) (17) 

- 1 ( e ) - 1 ( e )<^(e) L ( £ ) _ 1 (e^O) 1 (e) 



which are quantities that are free from singularities. For e = 0, they obviously include no contributions from interaction 
kernels K, and K» v while for e = 1 we recover Q {1) = G' 1 , Qfa = and = r^ oPoLE = ~ T^ OLE . With 
these definitions, one can easily verify that the interaction kernels can simply be written as 



%K = G (1) - G (0) (18) 
(1) y (o) 



iK^g^-g* (19) 



uc* = gft - gft (20) 



We will come back to the definitions above in section IV A since they obey Ward-Takahashi identities in a unified way 



due to their lack of poles; moreover, a derivation of Ward-Takahashi identities for the n-photon interaction kernels 
is straightforward with the help of eqs. (p8|)— (|2C|) . 
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B. Bethe— Salpeter Amplitudes 

With the systematic expansion in eq. (R), it is also possible to derive an "extended" Bethe-Salpeter equation 

X A = -iG A K A X A (21) 

for a bound state in an external field; here, we expand = Go,K,x- For the first order result, one can easily show 
that 

X 7 = GT^x (22) 

(see also [Q ) . With a little algebra, we find the second order result (see appendix |b|) : 

X 77 = GT 77 x . (23) 

Similar relations hold for the adjoint amplitudes if we start with the "extended" adjoint Bethe-Salpeter equation 
x A = _ i% A R A G A. 

X 7 = XT 7 G , (24) 

X 77 = %r 77 G . (25) 

For the sake of completeness, we finally want to give these results for distinguishable photons: 

X»(P',P)=G(P , )T»(P , ,P)x(P) , (26) 
X^{P', P) = G{P')T» V (P', P)x(P) (27) 

in first and second order, respectively, and analogously for the adjoint amplitudes x M (P',P) and x pl/ (P',P). In this 
sense, the expanded (adjoint) Bethe-Salpeter amplitudes x^ an d X M " (x M an d X M ") are n0 "new" quantities obeying 
separately relations like n-th order (adjoint) Bethe-Salpeter equations but they can be formulated in terms of the full 
Green's function, the full n-photon vertex and a "genuine" or "free" (adjoint) Bethe Salpeter amplitude without any 
coupling to the electromagnetic field. 

Let us stress that eqs. ( pq ) and ( p6| ) represent well defined relations since the Green's function G remains finite for 
q = P' — P 0; we will come back to this point later. 



IV. WARD-TAKAHASHI IDENTITIES 



Ward-Takahashi identities (see [|14 15 1) relate amplitudes of electromagnetic processes of order 0(e n+1 ) with those of 
order 0(e n ). In the following, we will derive constraints for the one-photon and two-photon irreducible interaction 
kernels from these relations as well as identities for the Bethe-Salpeter amplitudes in first and second order of the 
electromagnetic coupling. 

Since we aim at a general description of bound states composed of an unspecified number of ./V constituents, we will 
again in general suppress all relative momenta and indicate only the dependence of th e tot al four-momenta P and 



P'\ explicit examples of Ward-Takahashi identities for a qq system are given in section [V C 



A. Identities for the Amputated Green's Functions 



In p 

HP 



, a general form of the celebrated Ward-Takahashi identity is given (see also || and the seminal publications 

): 

q^{P',P)=e(M{P' -q,P)-M(P',P + q)) . (28) 



Here, M M denotes the amplitude for the process with n + 1 photons and M is the amplitude for the process with 
n photons; the indices for these n other photons as well as the dependency of M M and M of the internal relative 
momenta are suppressed in this notation. In fig. (|il|), this formula is depicted diagrammatically. 
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The amplitudes for processes with zero, one and two photons involved are given by 



M(P',P) 
M»(P',P) 
M^(P',P) 

respectively. The vertices (or amputated Green's functions) T' 1 and are given in eqs. (|12|) and (|13|). Further more 



G(P') T(P',P) G{P) 
= He) G(P') T»(P',P) G(P) , (29) 
= {-ief G(P') T^(P',P) G(P) 



we have defined r(J", P) := G 1 (P) ■ (2ir) 6 i (P' — P); for the charge factors in the equations above, see section VIA 



For the electromagnetic current in the order 0(e) with P' = P + q, the Ward-Takahashi identity reads explicitly 

g M M"(P + q,P) = e(M(P, P) - M(P + q,P + q)) 
<=> ~zq^(P + q,P) = G- 1 (P + q)-G- 1 (P) . (30) 

In the case of Compton scattering, four-momentum conservation requires P + q\ = P' + q 2 ] the Ward-Takahashi 
identity for this process in order 0(e 2 ) is therefore 

q^M^(P + qa, P) = e(M»{P - q 2 , P) - M V {P + qi -q 2 ,P + <fc)) 

<S=^ -iq lfi T^(P + 92 , P) = G-\P + q 2 )G(P - q 2 )T»(P - q 2 , P) 

-T%P + q 1 -q 2 ,P + q 1 )G{P + q 1 )G- 1 {P) . (31) 

It may seem surprising that not a simple difference of one-photon vertices occurs on the right-hand side of the last 
equation; however, we stress (as it will turn out in the next sections) that this combination of Green's functions 
and their inverse being attached to T v is crucial for a correct description of the Compton scattering process. A 
similar Ward-Takahasi identity in the order 0(e 2 ) has also been used in ref. [^2j in order to derive the lowest order 
contributions (i.e. at order 0(1) with respect to the photon momenta) for Compton scattering off spin-0 particles. 
The formulae in eqs. (^) and ( |3~l| ) provide general constraints for the explicit one-photon and two-photon vertices 
and since a violation of these Ward-Takahashi identities leads to electromagnetic observables that are not 
gauge invariant; we will show this in section 

Before we now insert the explicit expressions for the amputated Green's functions given in eqs. ( |l2] ) and (|l3|), we 
briefly recall the well-known identities for the free one-fermion propagator 

-i(P>-P)^ = Sp 1 (P')-Sp 1 (P) (32) 

and 

r = i-£rS?(P) (33) 

that — because Go includes only free propagators — immediately lead to the following relations for the lowest order 
vertices Tq and T^ v : 

- %qJ%{P + q,P) = Go^P + 9)- Go\P) (34) 
and - iqnJ% v (P + qi - q 2 ,P) = G Q \P + q x - q 2 )G (P - q 2 )K(P - q 2 , P) 

-r»(P + q 1 -q 2 ,P + q 1 )G (P + qi)G Q 1 (P) . (35) 
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We stress that due to the definitions in eqs. (£9|), the Ward-Takahashi identities in this section do not include any 
charge factors in physical units. However, these charge factors e and e 2 , respectively, can be re-introduced by a simple 



re-de finition of the n-photon irreducible vertices (see section VI A ) ; see also the remarks on charge factors in section 



IV C 



Let us finally make some remarks on the apparently different form of the Ward-Takahashi identities of order 0(e) 
and 0(e 2 ), respectively. The difference between the one-photon and the two-photon irreducible vertex is that 
in contrast to T M includes pole terms. To suppress the contributions from rp^ LE , we have introduced r> jn LE = 
r Miy — TpQ LE in eq. (|l4|). A generalization with and without inclusion of interaction kernels is given in eq. (|17]); it is 
straightforward to show that these quantities satisfy the following simple identities (see also fiQ): 

- iq^gfc (P + q,P)= G [t) (P + q)- e (e) (P) (36) 
znd -i qi ^ ) (P + q 1 -q 2 ,P) = g{ e) (P-q 2 ,P)-g{ e) (P',P' + q 2 ) . (37) 

Both for e = 1 and e = 0, they can be derived analogously to the Ward-Takahashi identities given above. We 
should note in general that all identities in order 0(e 2 ) can be re-formulated for the contraction with respect of the 
four-momentum q 2l , of the second (outgoing) photon by using crossing symmetry (i.e. q\ «-» —q 2 , /i <-> v). 



B. Identities for the Interaction Kernels 

As we will see in section we have to demand that the Ward-Takahashi identites in eqs. © and © hold in order to 
guarantee strict gauge invariance. From this principle, we can derive analogues identities for the irreducible interaction 
kernels. Let us start by inserting the explicit form of T^(P' , P) into the left-hand side of the Ward-Takahashi identity 
of order 0(e): 

- iq^{P + q,P) = ~iq^(P + q,P)~ iK^{P + q, P)) 

= G- 1 (P + q)-G 1 (P)-q^K^P + q,P) . 

We see from eq. ( |30| ) that the difference of two inverse full Green's functions should appear; from this constraint and 
due to G _1 = Gq + iK, we can read off a Ward-Takahashi identity for the one-photon irreducible interaction kernel: 

-i qfi K»(P + q,P) = K(P)-K(P + q) . (38) 

This result is well-known for a long time and has been derived e.g. in refs. [^,0. However, the same procedure can 
also be applied to the divergence of the two-photon vertex T^(P' , P) with P = P + qi — q 2 and yields 

- iq llM T^(P + qi - q 2 ,P) = -iq lfl T^(P + q x - q 2 ,P) iq lfl K^{P + Ql - q 2 ,P)) 

+T"{P + qi -q 2 ,P + qi ) G(P + qi)(- iqi^(P + q u P)) 
-T V {P + qi -q 2l P + q X )G {P + <&)(- iqi^(P + q u P)) 
+ ( - iq^T^P + qi -q 2 ,P- q 2 )) G(P - q 2 )T» [P - q 2 ,P) 
-( - iq lfi T%(P + qi - q 2 ,P - q 2 ))G Q (P - q 2 )T»{P - q 2 , P) 
= G- l {P + qi -q 2 )G{P- q2 )T v {P-q 2> P) 
-T V {P + qi - q2 ,p + qi )G(P + qi )G-\P) 
+r»(P-q 2 ,P)-r»(P-q 2 ,P) 

-r»(p + qi - q2 ,p + Ql ) + T V {P + Ql -q 2 ,P + qi ) 

- qifi K^(P + qi -q 2 ,P) 

where we have used eqs. (fjOl), (|34|) and (35). By comparison with eq. ([H]) we observe that the correct terms for the 
Ward-Takahashi identity in 0(tr) are already present and with Tq — T v = iK u find the relation for the two-photon 
irreducible interaction kernel: 

-iq^K^(P + q 1 -q 2 ,P)=K v (P-q 2 ,P)-K v (P + q 1 -q 2 ,P + q 1 ) . (39) 



Note that we have assumed that K v obeys the relation in eq. (|3§|) since we have used the full Ward-Takahashi 
identity in the order 0(e). Let us remark that the identities in eqs. ( |38| ) and (^) can be derived in an even more 
transparent way by using eqs. (^8|) - (|20|) together with the Ward-Takahashi identities for the non-singular quantities 
0^ and Q?K in eqs. (|(|) and 
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C. Example: Explicit Identities for qq Mesons 



In eq. (f28[), we have omitted the dependence on the relative momenta in the amplitudes M M (P',P) and M(P', P); 
in general, a sum over all photon couplings to the internal fermion lines and therefore a dependence of all internal 
momenta should be written out. We will demonstrate this for a mesonic system explicitly, i.e. for a bound state 
composed of a quark q with momentum p\ and charge e\ and an anti-quark with momentum — p 2 and charge — e 2 . 
The two-body Ward-Takahashi identity in this case reads (see jl , (|] and references in || ) 



iq^M t *(p' 1 ,-p' 2 ;p 1 , -p 2 ) = e 1 (M(p' 1 - q, -p' 2 ;p 1 ,-p 2 ) - M(pi, -p' 2 ;Pi + <Z, -P2)) 
- e 2 (M(p[, -p' 2 + q;px, -p 2 ) - M(p[, -p' 2 ;pi, ~Pi - q)) 



(40) 



note the details in the coupling of the photon to the anti-quark. For the sake of simplicity, we assume that the 
operation of the the charge operator Q/ is trivial in flavour space and only produces a flavour dependent factor (e.g. 
e.i = — l/3e or = +2/ 3e), i .e. that Qj is diagonal in the fi = u,d,s flavour basis and acts like Q/|/i) = ei\fi); see 
also eq. (|7(]) in section VIA . We shall assume here that the interaction will not be mediated by charged exchange 
particles (i.e. [K, Qf] =0). Furthermore, we define the total b oun d state momentum by P = p\ +p 2 and the relative 
momentum by p = r\ 2 p\ — r\\p 2 (with 771 + 572 = 1 ) , see also eq. ( |A6| ) . With P' = P + q, the two-body Ward-Takahashi 
identity in this new coordinates can be written as 



i^M^PV^p) = ei(M(P' -q,p'- mq;P,p) - M (P', P '; P + q,p + V2 q)) 
- e 2 (M (P 1 ~q,p' + m q; P, p) - M (P',p'; P + q,p- Tfc?)) 



(41) 



Now let us check the free Green's function proportional to the free amplitude Mg 
written out in terms of one-particles propagators: 



-ieGp = — ieGoT^Go, explicitly 



iMg(P\p';P,p) = eG (P', P ') [T% tl (P',p' '; P,p) + T^ 2 (P' \p' '■ P,p)] G (P P ) 

— 1 



(42) 



My +p')®s f (-^+ v ') 



ei ■ 



e 2 ■ 



P 



pi 



S^(-^+p)\-(2nfS 4 (P' 



F ( Y+ P') 



1 y 



P - q)S 4 (p' - p- i] 2 q) 

(2n) s 6 4 (p' -p-q)S 4 ( P ' -p + mq) 



Sf( j + P) ® S F (-- 



p) 



where Tq couples the photon to each fermion line of the free bound state propagator Go- Using the one-particle 
Ward-Takahashi identity in eq. (p3), we thus find 



-iq„M^{P',p';P,p) 
— —ie\ 



(43) 



S F (-+p)®S F (-- 



-ie 2 



S F (-+p)®S F (-- 



■p) 
■p) 



Sf(^+p')®S f (-^+p) 

s F (L +p) ®s F (~+p') 



(2^f5\P'-P-q)5 A (p l -p-i l2 q) 
(2vf5 4 (P' -P-q)6 4 (p> -p + mq) 



We now define §j :— ej/e and write down the related expressions for the amputated one-photon irreducible vertices 
eT» = iG-^M^G- 1 and e r^ = iG^ 1 M^Gq 1 : 



iq^V (P',p'; P p) = St (G- 1 (P', p'; P + q,p + r) 2 q) - G^P' - q, p' - r) 2 q; P, p)) 
-e 2 (G- 1 (P , ,p';P + g,p-77 1( 7)-G- 1 (P , -q,j/ + m ;P,p)) 



and 



iq ll T%(P',p';P,p) 



e^G^iP^p'; P + q, p + mq) 
~e 2 (G^(P\p'-P + q,p-i liq ) 



G \P' 

here, the inverse free propagator Go (depending on primed and non-primed coordinates) is defined by 



q,p 
q,p' 



V2q;P,p)) 
mq; p ,p)) 



(44) 



(45) 
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G^(P',p';P,p) 



(2tt) 8 <5 4 (P' -P)5 A (p' -p) . (46) 



If we recall G 1 = G 1 + iK as well as the general relation = Tq — iK^ for the one-photon irreducible vertex, 
then we easily find with P' = P + q 



iq^K" (P',p'; P, p) = h [K(P' -q,p'- V2 q; P,p) - K(P' ,p' ; P + q,p + m )) (47) 
- e 2 (K(P' -q,p' + mq; Pp) - K(P',p'; P + q,p- rnq)) 



This is the full version of the constraint in eq. (pq), formulated explicitly for a qq system without omitting the 
relative momenta. As we have pointed out in the last subsection, the divergence of K 11 vanishes for q — > 0. We now 
see in detail that this also holds if the interaction kernel K is independent of the total four-momentum P and is 
simultaneously of local type; then the terms on the right-hand side of eq. ( f47j ) mutually cancel if K ~ K{p' — p) is a 
convolution-like kernel. 

In analogy to this derivation, we find the explicit formulation of the Ward-Takahashi identity for the two-photon 
irreducible interaction kernel of a qq system with P + q\ = P' + q 2 : 



iq^K^ (P' ,p';P,p)=e 1 [K v (P' - qi ,p' - rj 2qi ; P,p) - K V (P' ,p'; P + q uP + mqi )j (48) 
- e 2 (K V {P' -q 1 ,p' + r ll q 1 ;P,p)-K u (P',p';P + q 1 ,p- mqi ) 



Note that for a deuteron described as a pn bound state (i.e. with e\ — e p — 1 and e 2 — e„ = 0), the r.h.s. of this 
equation reduces to the difference in the first line. The resulting identity (up to a factor ie and in a slightly different 
notation) for the relative momentum parameter choice rj\ = rj 2 = \ is completely equivalent to eq. (4.13) in ref. 
fllf ; in this publication, the authors have found this constraint for the two-photon interaction current by demanding 
gauge invariance of the Compton scattering amplitude in a general Bethe-Salpeter framework. In fact, this represents 
a completely different approach to this issue compared to the one shown here. While we find the Ward-Takahashi 
identities for the n-photon irreducible interaction kernels in a rather general framework starting from the "gauging 
of equations method" applied to the Green's function of a composite system, the authors of ref. Jll| begin with the 
description of a special bound state (i.e. the deuteron system); however, they keep a certain universality since they 
do not specify the interaction kernel in the Bethe-Salpeter equation. Then they consider Compton scattering off this 
system, include the lowest order contribution (~ Tp"), re-scattering terms via T matrix diagrams and one-photon 
interaction kernels K^, and study the gauge invariance of the resulting tensor. They find that —iqi^T^ — only 
holds if an additional two-photon irreducible term is included; this condition is then recasted in the condition quoted 
in eq. (j4|). 

The direct derivation of the identity in eq. ( pl8| ) is lengthy since for instance the (explicitly crossing-symmetric) 
subtraction terms in eq. ( |l3| ) for a qq system read explicitly 



(49) 



r£(p', p' + q 2 )G (p' + q 2 )r£(p + qil p) + rft(p', p' - qi )G (P' - gi)r£(p - q2 , p 

= e\ Tl x (P',p' + V2 q- P 1 + q 2 ,p" + V2qi )Go(P' + q 2 ,p" + r^q^P + q U p'" + V2QlK A (P + q x ,jf" + Tfeft; P,p) 
+ el ro j2 (P',p' ~ViT,P' + Q2,p" -rjxqi)G (P' + q 2> p" - T]iqi\ P + qi,p"' - rjxqi)T^ 2 (P + qx,p'" -rjxqi;P,p) 
- e 1 e 2 r- 2 (P',p'+ Q ;P' + g 2 ,p'' + 7 72 g 1 )Go(P' + g 2 ,p'' + ^gi;P + gi,p''' + ^gi)r(; i (P + ( Zi,p' , ' + ? 729i;P,rt 

i g 2 ^o,i(p\p - Q ;P' + Q2,p" - i]in)Go(P' + q 2 ,p" - mqi, p + qi,p"' - vm)roAP + h,p" - mn\P,p) 



qi <-> -?2 

where we have defined q := qi — q 2 and Q := r\ 2 q\ +r/iq 2 ; the definition of Tq i can be read off eq. (|42|). Note that we 
still neither indicate the integrations J cftp" and J cftp'" nor the various indices; see eqs. (|l|) and ((^) for comparison. 
Due to the complexity and the large number of terms occuring in the explicit qq formulation of the two-photon 
irreducible vertex , we feel that it might be more instructive to keep on omitting the relative momenta and the 
terms related to internal charges in the rest of this contribution since it will clarify the structure of the resulting 
expressions. We stress however that from contributions like those quoted above it is in fact straightforward to derive 
the explicit Ward-Takahashi identities for a qq system in the order 0(e 2 ). 
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We have chosen the explicit example in eq. (E^) not only to justify the suppression of relative momenta in this 
paper but also to illustrate the role of the subtraction terms — T^GqTq — TqGqT^. Note that if we write out the 
crossing symmetric terms in this example by exchanging q% <-> — q% and /i <-> v and by explicitly using the one-particle 
propagator formulation as in eq. ([l2|), we find that the ef parts are correctly produced but that a doubling of terms 
proportional to e\&2 occurs; these terms are essentially proportional to 7^ ® 7" and 7" ® 7^, respectively. This 
is precisely what happens in the lowest-order parts of +r M GT !/ + T^GT^ — only with an opposite sign. By this 
mechanism, it is guaranteed that the only lowest-order terms surviving in the complete two-photon irreducible vertex 
r M " are collected in Tq V without any double-counting. 

D. Identities for the Bethe— Salpeter Amplitudes 

Now that we have derived the Ward-Takahashi identities for the n-photon vertices and irreducible interaction kernels, 
we are able to give similar relations for the "extended" Bethe-Salpeter amplitudes of orders 0(e) and 0(e 2 ). 
If we take the divergence of eq. (^) , then we find immediately with the Ward-Takahashi identity for the one-photon 
vertex the relation 

-iq^(P + q,P) = X(P) ; (50) 

here, we have used the Bethe-Salpeter equation G~ 1 (P)x(P) = 0. Analogously, the relation for the second order 
Bethe-Salpeter amplitude reads 

-i qi ^(P + qi-q2,P) = x' y (P-q2,P) ■ (51) 
As we have already stated in section III B| , the relations above are well defined for finite photon momenta. 



E. Differential Identities 

We will now consider differential forms of the relations in the preceeding subsections. The limit q,q\,q2 — > yields 
P' — > P; therefore we suppress all momentum dependencies in the following. 

With := d/dPfj, or equivalently assuming infinitcsimally small photon momenta, we find for the extended Green's 
functions 

d^G = tG" (52) 
and = iC" . (53) 

From eqs. (|3^) and (|39|), we obtain the following relations for the one-photon and two-photon irreducible interaction 
kernels: 

d^K = iK^ (54) 
and d v K» = iK» v . (55) 

In order 0(e), the situation for the amputated Green's function is simple. However, in the second order relation in 



eq. (31) no simple difference with shifted arguments occurs on the right-hand side. To derive a Ward identity for 



I 1 ^, we therefore start with eq. flSJ) and recall that G^ = GT^G: 

= (d u G)r>*G + G{d"T^)G + Gr»(d v G) 
= {iG u )T»G + G{d v T»)G + GT^(iG l/ ) 
= G(iT u GT» + iT^GT" + d v T^)G 

= iGT^G . 

For the one-photon and two-photon vertices, the corresponding relations with Y := G~ x are therefore 

= -iT^ 1 (56) 
and d v T» = -i(T v GT>* + T^GT V - T^) . (57) 
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Again we find that the differential two-photon relation for the vertex is not a simple analogue to the one-photon case; 
however, both orders of the electromagnetic coupling have in common that exclusively non-singular terms appear since 
we can write ST* 1 = ^r^c-PoLE ^ or ^ ne l &s t identity, see eq. We will see in section fvi| that the non-singularity 

of the last equation will be crucial for the correct low-energy limit of the Compton scattering tensor; note that the 
sign in eq. ( CT ) is essential for the related cancellation in the limit qi — ■> 0. 

For the sake of completeness, we also quote the analoguous results with Tq := Gq 1 for the free vertices: 



and 



9T = 
0Tg = 



-*rg 



-i(r G Q y$ + r%GoTZ-r$') 



(58) 
(59) 



We will finish this section by studying differential identities for \i X^ and X^ v comparable to those above. Let us 
anticipate that the following relations basing on eqs. (^) and (^) are derived in a quite formal way, i.e. regardless 
of the poles that will appear in the Green's function G{P) for P — > P (P 2 — M 2 ) or equivalently for vanishing photon 
momenta 



To find the derivative of the Bethe-Salpeter amplitudes, we cannot start from the relation in eq. (50) since there no 
difference appears. However, we can exploit the Bethe-Salpeter equation and find 



d^~iG K x ) = -i((d»G )Kx + G (d»K) x + G Q K{d» X )) 
(1 + iG K)-\ Go (rg GoKx +K» X ) 



(60) 



-GG- 



= iG{T» ~ iK») X = iGT^x ■ 

The second derivative of x can be calculated analogously by using eq. 
occur; we find with eqs. (E6I) and (E7j) the following relations: 



and d v X tl = ix^ = iGJT'x 
Similar equations hold for the adjoint amplitudes: 



and 



d»x = ir =ix^G 
d v x lt = ix^ v = iX^ u G 



57|) such that characteristic cancellations 

(61) 
(62) 

(63) 
(64) 



Note that the bound state Bethe-Salpeter amplitudes are defined only for on-shell momenta P 2 = M 2 whereas the 
Green's function G(P) generically includes off-shell contributions. Therefore the above-mentioned equations actually 
have to be considered with care: any derivative of x would infinitesimally shift the Bethe-Salpeter amplitude away 
from mass-shell in a region where x becomes ill-defined. However, we can formally introduce a quantity X (P) for 
off-shell momenta P 2 ^ M 2 analogously to the definition of x m cc l- ©! then X(P) is not a solution of the Bethe- 
Salpeter equation but we can properly define the operation d^X. Additional terms due to G~ X X ^ will appear in 
eq. ( |60|) that vanish in the limit P 2 — > M 2 in which we find the final relations in eqs. (|6l|)-(|64|). 
We shall remark here that the Green's function G that enters in eqs. (|6l|)-(p4) exhibits a pole in the on-shell case 
and (for states with total momentum P and mass M) can be written (see e.g. [131) as 



G(P) 



X(P) ° X(P) 
P 2 - M 2 



terms regular 
for P 2 -> M 2 



(65) 



Note that the Bethe-Salpeter amplitudes are on-shell quantities and therefore depend on P with P 2 = M 2 . The 
physical implication of the dominant term in the Green's function G(P) for the derivative of the Bethe-Salpeter 
amplitude and the related first order expanded amplitude x^ m the limit P 2 — > M 2 is depicted diagrammatically in 
fig. (H|); an analoguous diagram can be sketched for the second order Bethe-Salpeter amplitude x^ ' ■ 



d =] 
dP u 



x r 




X 




P 2 -M 2 



, regular 
x r ~r terms 



FIG. 12. The Bethe-Salpeter amplitude expanded in first order of the electromagnetic field according to eq. (|si]). 
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We stress again that these differential relations for the Bethe-Salpeter amplitudes \ an d X are derived with quite 
formal arguments and without discussing the apparent poles that enter the equations above. However, they allow for 
an alternative derivation of low-energy limit for the Compton scattering process as we shall show in section VI C . 



V. GAUGE INVARIANCE 



It is an important question whether the one-photon and two-photon vertices defined in the eqs. (|12|) and (13) satisfy 
gauge invariance if we evaluate them between Bethe-Salpeter amplitudes. To formulate this task more precisely, let 
us first define the electromagnetic current via 

J" := (P' |j"(0)| P) = -x(P')T»(P', P)x(P) (66) 

and the Compton scattering tensor via 

I*" : = i [ d 4 x(P' \Tj»(x)f(0)\ P)e- iqix = ix(P')r^(P', P)x(P) ■ (67) 



Here, T denotes the time-ordered product. Gauge invariance implies that the divergences of J M and T^ v vanish, i.e. 
that in momentum space 

- ig„J" = and - iq lfl T^ = -iT^q 2v = 

holds. In the subsequent paragraphs, we will address this question for one-photon and two-photon amputated Green's 
functions in various approximations. Let us furthermore stress that in the following two sections the momenta P and 
P' of the incoming and outgoing bound state are supposed to be on-shell (i.e. P 2 = P' 2 = M 2 ) since otherwise the 
Bethe-Salpeter equations could not be applied. 



A. Electromagnetic Current 

With the Ward-Takahashi identity for the one-photon vertex in eq. ([30j), it is easy to show that the divergence of 
the electromagnetic current vanishes (here, energy-momentum conservation requires P' = P + q): 

- iq„J" = -X(P')(- k^(P', P)) X (P) 

= - X (P')(G- 1 (P')-G- 1 (P))x(P) = . 

Obviously it is trivial to show the gauge invariance of the electromagnetic current with the help of the Bethe-Salpeter 
equations G~ 1 (P)x(P) = xiP')G~ 1 (P l ) = 0. However, this requires that the one-photon vertex obeys the Ward- 
Takahashi identity in eq. (pBG) which is only true if the one-photon irreducible interaction kernel satisfies the 
relation in eq. (|3^). To see what happens if we neglect contributions of this one-photon kernel, i.e. if we only 
consider the lowest order vertex Tq, we again exploit the Bethe-Salpeter equations Gq 1 (P)x(P) = —iK(P)x{P) and 
X {P')G^{P') = -ix(P')K(P') and find 

- i%J$ = -X(P')(- iq^(P',P)) X (P) (68) 
= - X (P')(Go 1 (P')-G^(P))x(P) 
= ix(P')(K(P')-K(P)) X (P) ■ 

This expression only vanishes if the interaction kernel does not depend on the bound states' four-momenta. Further- 
more it turns out that the kernel has to be of local type — i.e. K(P; {kj}, {kj}) — K({kj — k'j}) — if we explicitly 
write out the dependence on the relative momenta. This is true for kernels in ladder approximation but not necessarily 
for instantaneously approximated kernels (as e.g. used in [ p~8|Jl9| ) . 

Note that, independent of the kernel type, we can state that —iq^J^ vanishes and therefore gauge invariance holds 
for <? = P' - P ^ 0. 
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B. Compton Scattering Tensor 



As for the electromagnetic current, it is again easy to show that the divergence of the Compton scattering tensor 
vanishes if we take into account all contributions in T^"(P' , P). Then the Ward-Takahashi identity for the two-photon 
vertex in eq. (31) holds if the included n-photon irreducible interaction kernels obey eqs. ( j38|) and (|39|). With the 
Bcthc-Salpcter equations, we therefore find for P 1 = P + qi — q 2 

- iq Xtl T» v = ix(P')( - iqi^(P>, P)) X (P) 

= ix(P')(G-\P')G(P - q 2 )T v {P - q 2 ,P) - r"(P', P + q x )G{P + (Zl )G- 1 (P)) X (P) = . 

As an approximation, we now want to neglect the contributions from one-photon and two-photon irreducible inter- 
action kernels. Then the subtraction term in cq. ([l3]) cancels the lowest order part in r 7 GT 7 — see eq. ([IT]) — and 
we start with a modified two-photon vertex reading explicitly 

f^(P', P) = T^(P', P) - iT^P', P' + q 2 )T{P' + q 2 )T%(P + q x ,P) 

-iT^P',P'-q x )T(P -q 2 )T%(P-q 2 ,P) , (69) 

where the T matrix satisfies eq. (||). The divergence of this vertex can easily be derived from the Ward-Takahashi 
identities for the free one-photon and two-photon vertices in eqs. (|34| ) and (p5|) : 

- iq^f ^(p', p) = g \p + qi - q2 )G (p - <&)rg(P - q 2 , P) 

~ T%(P + qi -q 2l P + q x )G (P + qi )Go X (P) 

- iI*(P + qi -q 2 ,p + qi )T(P + gi)(Go X (P + q x ) - G J (P) 

- i(G \P + qi - q 2 ) - G„\P - q 2 ))T(P - q 2 )r»(P - q 2 , P) . 

Sandwiching this expression between a Bethe-Salpeter amplitude and its adjoint, using Gq x — an( i X.Gq 1 = 

— ixK, respectively, and inserting the T matrix equation — see fig. (||) — , we obtain the divergence of the Compton 
scattering tensor as 



iq Xfi T^ = ix(P') [ - iqx^{P\ P)Jx(P) 

= ix(P') (i-iK(P'))G (P - q 2 ) - i{-iK{P>))T{P - q 2 ) 

+ iG^(P - q 2 )T(P - q 2 ))T v (P - q 2 , P) X (P) 
-*X(P') W, P + q x ) (Gq(P + qx){-iK{P)) - iT{P + q x )(-iK(P)) 

+ iT(P + qi )G Q \P + q 1 ))x(P) 
= x{P')(k{P')G {P - q 2 ) - iK(P')T{P - q2 ) - K(P - q 2 )G Q (P - q 2 ) 



+ lK{P - q 2 )T(P - q 2 ))T»{P - q 2 , P) X (P) 
-x(P')r£(P', P + qi )(G (P + qi)K(P) - iT{P + qi)K(P) - G (P + qi)K(P + q x ) 

+ iT{P + q l )K{P + q 1 ))x{P) 
= x(P')(K(P') - K(P - g 2 )) (Gq(P - q 2 ) - iT(P - q 2 ))r»(P - q 2 ,P) X (P) 
+ X (P')n(P', P + qi )(G (P + qi) - iT{P + qi )\ (k(P + qi ) - K(P))x(P) ■ 



With G = G - IT and defining x' 1 := GT^x (where V = r£) in analogy to eqs. (|22|) and (|23|), we can therefore 
write the divergence of Compton scattering tensor without inclusion of n-photon irreducible interaction kernels as 



iqi^ = x(P') (K(P') - K(P' - qi )) r (P - q 2 , P) 
+ f(P', P' + (ft) (k(P + q x ) - K(P))x(P) 
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This expression vanishes — as in the similar case of the lowest-order electromagnetic current, see eq. (|6S|) — in 
the limit qi,q 2 — > (P' — > P). The gauge invariance condition —iq\^T llv = also holds if we use an interaction 
kernel K(P) that does not depend on the total four-momenta of the bound states. As it has been pointed out in the 
preceeding subsection for the lowest order current Jq, an inclusion of all relative momenta in the calculation shows 
that the kernel must also be local in order to satisfy gauge invariance in this case. This result has also been found in 
p| where an explicit calculation for deuteron Compton scattering has been presented. We agree completely with the 
authors who stress that this observation emphasizes the necessity of the inclusion of n-photon irreducible interaction 
kernels if the constituents of the composite system interact via non-local potentials. 

The last point to study is the question of gauge invariance for the lowest order two-photon vertex Tq U . If we recall 
the free Ward-Takahashi identity for this vertex in eq. ( |35| ) , we find immediately 

- = ix(P'){ - iqx^{P',P)) X {P) 

= zx(P')(G ( 7 1 (P')G , o(P - q 2 )T v (P - q 2 ,P) - T»(P', P> + qi )G (P + qi )G^ (P)) X (P) ■ 

Without applying the Bethe-Salpeter equation, we see that this lowest order expression only vanishes for qi, q 2 — > 0. 
It is clear that in this limit the re-scattering terms due to the T matrix will not contribute; therefore gauge invariance 
in the very-low-energy regime will be approximately restored even without explicit inclusion of intermediate states. 



VI. LOW— ENERGY LIMITS 



In this section, we will study the low-energy behaviour of the one-photon expression xT^X an d the two-photon 
expression xT^X- We will restrict our discussions to bound states with total angular momentum J — so that we 
will not have to consider technical subtleties originating from non-zero spins; however, our statements remain valid 
also for spin-averaged amplitudes describing Compton scattering e.g. off a proton or a neutron with total spin i. 
In the first subsection, we will give a concise formalism how to include the correct charge factors. To clearify this 
point, we will pin down the lowest order results for the electromagnetic current and the Compton scattering tensor 
in explicit expressions for the bound state being a (pseudoscalar or scalar) qq meson as an example. 



A. Charge Factors 




If we now study the low-energy limits of the electromagnetic current and the Compton scattering tensor, a caveat 
is in order at this point. Up to now, we have considered the photon coupling in Dirac space only, which was of the 
(correct) vector type 7 M ; we thus neglected the (in the general case iso-spin dependent) charge operator Q/ in flavour 
space that reads e.g. for hadronic bound states 

° \ 

(70) 

-U 

for constituents being u, d and s quarks (see also for a different formulation of this issue). This observation 
does not spoil our results on Ward-Takahashi identities and gauge invariance but it turns out in this section that it 
will become important for the low-energy limits. Instead of the operator in eq. ( f70| ) , we simply get Q f = 1 if we 
introduce our vector coupling in Dirac space by a derivative of the one-particle propagators. However, these fermion 
propagators Sp(Pi) are functions of the momentum of the i-th constituent which depends in a simple way of the total 
momentum P (see appendix ^). This dependence is given in eq. ( |A5| ) and yields 

Q 

Sp(pj) = ~ir) 3 S F {pj)^S F {pj) . 



dP» 



As described in the appendix, the rjj's are parameters that set a special choice of relative momenta; we can fix them 
by Vj = e j/Q w ith Q = J2j 6j being the total charge and ej being the charge of the j-th constituent. 
We now have to specify what we mean in detail with the lowest order one-photon vertex Tq for a composite system 
of N constituents. Recalling Gq 1 —: To where the total momentum is unchanged by Go (i.e. P = ^2jPj = ^jPj)^ 



we state with eq. ( |58| ) that 
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i—Go 1 (P,{k! i },{k i })=i—Go 1 ({ P '},{Pi}) 



d _d_ 

^({^■},te}) = E r o,i({^}'te}) 



ap M u v dp^ 



= m-l li ® Sp 1 ( P2 ) ® S^ X (P3) ® •• • ® ^(pjv) • <5i({^- -p 3 -}) 
+ 772 • S^fai) <g> S^ips) ® ... ® S^^pjv) • * a ({pJ -Pj}) 
+ »73 • SfHpi) ® ^(Pa) ® y ® ••• ® S^CPJv) ' *»({p$ - Pi}) 

+ ?7w • S£ X (pi) ® Sp 1 { P2 ) ® S^ x (p 3 ) ® • ■ ■ ® 7 M ' M{Pj ~ Pi» 
with £ = 1,2,..., iV — 1 and j,k = 1,2, ... ,N and the short-hand notation 

SMi ~ Pj}) ■■= (2tt) 4 ^ ■ $Vi - P1 )6\ P ' 2 - p 2 ) . . . <5 4 (p' fc _i - p fc -i)<J Vfc+i - p fc+1 ) . . . 5\p' N - PN ) 
The coordinate choice ?7j = ej/Q induces a re-definition of the free vertex in eq. (|7l]) by 

N „ JV 



J- n 



o=E r o, — ^ r o=^E r o, ( 71 ) 



such that now 

= ej • Sp X { Pl ) ® . . . ® S^fe-i) ® 7 M ® ^(Pj+i) ® • •• ® Sp X {p N ) ■ Sjdp'j - Pj }) 

holds. This procedure can also be applied to the full one-photon vertex where also contributions from the one- 
photon irreducible interaction kernel K M are included; in all terms that are considered the photon can be coupled to 
each fermion line in the corresponding diagram. Analogous re-definitions for the two-photon vertices, i.e. 

> J_p^ I/ (72) 

can also be given. We stress that this scheme is essentially equivalent to the introduction of a proper electromagnetic 
charge operator by hand followed by an evaluation of the contraction of all indices in flavour space. 
In the following subsections, we will give — together with the general results — also the explicit lowest order formulae 
for qq mesons to make clear that these re-definitions give the correct and complete matrix elements for electromagnetic 
form factors and Compton scattering. 



B. Current Normalization 



From Lorentz invariance and charge conjugation symmetry, one can deduce the matrix element of the electromagnetic 
current to be of the on-shell form J M = Q ■ (P + P')^f(q 2 ) where f(q 2 ) is the form factor with the normalization 
/(0) = 1 and Q is the total charge of the bound state. In the limit of vanishing photon momentum q = P' — P — > 0, 
we therefore find 

lim J* 1 = lim (P' | f (0) | P) = Q ■ 2P fl . 

q — >0 (J — >0 

Let us now consider the normalization condition of the Bethe-Salpeter amplitude describing a bound state of mo- 
mentum P and mass M (see e.g. fjjJ): 



X{P) {^p jG- \P) + iK{P)X ) | p3=M2 X (P) = i 2P» 



=G-i(P) 

With the Ward-Takahashi in eq. (|56|), recalling G _1 =: T and by comparison with the definition of the current in 
terms of Bethe-Salpeter amplitudes and the one-photon vertex — see eq. ( |66| ) — , we easily find the correct form 
factor normalization for q — > 0': 
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(73) 



The total charge Q on the right-hand side comes from the re-definition of the vertex in eq. (|7lj) acting such that a 
modified Ward-Takahashi identity d^G^ 1 = —i-^T^ 1 has to be applied. 

This low-energy limit is even correct if we neglect contributions from the one-photon irreducible interaction kernel . 
Due to the Ward identity d^K = iK M , the interaction kernel K then must be independent of the total momentum P 
(as we have also found for the gauge invariance of the lowest order current Jq). We thus can apply the Ward-Takahashi 
identity in eq. ( |58| ) for the free one-photon vertex: 

-X(PK(P,P) X (P) = Q-W ■ (74) 

As an example, we give the left-hand side of this last expression explicitly for a qq bound state and a non-vanishing 
photon momentum q =/= (note that here the anti-quark is defined with momentum —p 2 and charge — 62): 

+ e 2 J -0^x(P,P^Viq)(s F 1 (viP + P)®r)x(P,p) • 

The charge factors introduced in the preceeding subsection are obviously correct; the indices of xXoX therefore must 
only be contracted in Dirac space (or, to be precise, the contraction in flavour space will yield a factor 1 due to the 
unaltered charge operator Qy = 1) . 



C. Born— Thomson Limit 



For charged particles, the Compton scattering tensor includes not only structure dependent terms but also pole 
contributions in order O(cu ) with u> = (qi,q 2 ); then the Born part of in soft-photon approximation is known to 
be (see e.g. §§) 



rpllV 

Born 



Q 2 ■ /(^)/fe 2 ) 



(2P + qi y(2P> + q2 y , (2P 1 - qi y{2P-q 2 ) 



M 2 



— M 2 



(75) 



with the Mandelstam variables s — (P + qi) 2 = (P' + q 2 ) 2 and u = (P — q 2 ) 2 = (P' — qi) 2 and the one-photon 
form factors f(qf) defined as in the preceeding subsection. Note that T-q^ is explicitly gauge invariant due to the 
non-singular term proportional to g^ v \ the three terms in eq. ( f75| ) can be identified with the diagrams (a), (b) and 
(c) in fig. (|l3]). The factor Q 2 forces the amplitude to vanish for neutral bound states in the low-energy limit. 






(a) (b) (c) 

FIG. 13. Born terms for Compton scattering off a spin-0 particle. 



In the non-relativistic limit (i.e. for — + 0), only the Born terms in the full tensor T^ v survive. The amplitude 
T for Compton scattering off a bound state of mass M with spacelike photon polarization vectors e% = (0,£i) (and 
therefore qi ■ £j = and P ■ E{ — P' ■ e.- L — for real photons) then yields the celebrated result of Low, Gell-Mann and 
Goldberger (see |l||||2§): 



Um n T =-VK7 Um n £ ^ T ^ £2 » = ~% (^1 ' ^) ■ ( 76 ) 

— ^0 AM qi^O M 

Obviously, the pole terms in eq. J75| ) do not contribute to the Compton amplitude in the very-low energy limit; note 
furthermore that here and in the following we have used /(0) = 1 for real photons. 

Now we want to show that the full two-photon vertex in fact yields this limit for qi,q 2 — > {P 1 — > P). Let us 
therefore consider eq. ( |73| ) and make some remarks on the right-hand side of this current normalization. Formally 
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the factor 2P M originates from a differentiation of the free boson propagator Ap(P) — (P 2 — M 2 ) 1 and a subsequent 
amputation of this vertex in the limit q — ► (see e.g. 



J" := lim J^ 1 = A F X (P) 



Af 1 (P) = Q-^A F 1 (F) = Q.2F^ . (77) 



The result is an amputated (point-form) boson-photon vertex; if we once more take the derivative — Qd/dP u then 
we end up with a two-photon vertex of non-singular form: 

nToP OLE == ^o T nopo LE = -2 • = -Q 2 • w v jp-*?iP) = -Q 2 - W ; (78) 



see diagram (c) in fig. (p3|). ^noPole * s equivalent to the non-singular part of Ijjorn * n ^ ne nm it <Zij <jf ~^ 0, see eq. 
( f75| ) . It can be recovered in the Bethe-Salpeter formulation by taking the derivative of the left-hand side of eq. ( |73| ) 
and applying —Q,d/dP v only to the amputated (point-form) two-photon vertex: 

-x(P)(-Q^rT»(P,P))x(P) (79) 
= lim ix(P') (T»»(P', P) - T U (P', P' + q 2 )G{P' + q 2 )T»{P + q U P) 

- r"(P', P' - qi )G(P - q 2 )T v {P - <h,Pj)x(P) 
= ]imix(P , )^ 0P0LE (P',P)x(P) 



Here we have used that c?T M = -Q^^ p 0LE , see eq. ( p7\ ) and eq. (72) for the vertex' re-definition. In analogy to eq. 
(|l4|), we now split the Compton scattering tensor into two parts by defining = Tp^ hE + ^noPole wner e PpQ hE 
includes all pole terms. Comparison with eq. ( |67| ) then yields 

^ T NoPo LE - ^ ix(P'm PohE (P',P) X (P) = -Q 2 ■ 2«T (80) 

which is the correct low-energy limit for the non-singular terms in the Compton scattering tensor. In the following, 
we will only work in the limit qi — > (i.e. P' —* P), therefore skip the dependence on the total four-momentum P 

and use the abbreviation d" := 75%-. 

aP v 

Let us now study the pole contributions in the low-energy limit of T Ml/ which are necessary ingredients of the full 
tensor with respect to gauge invariance. First, we will investigate the second line of eq. (79) and find that the terms 
Y^GY" + Y V GY^ indeed reveal the correct pole structure. This is due to the fact that the full Green's function G(P) 
has a pole for P 2 — > M 2 (and therefore in the limit of vanishing photon momenta) with the residue — %x X, see eq. 
(|65|). Recasting the whole expression and using eq. (18(3) yields 

Jft, (* r "* ' jhf 2 ■ X ^ X + X ^ X ■ ' xT " X + T noPo LE ) = }™ ixF^X ■ (81) 

Since J* = -x{P')Y^(P' ', P)x(P) = Q{P + P'Yfiq 2 ) with q = P' - P holds for the electromagnetic current, we 
recover the full Born term T^ RN of eq. ([75|) in the limit qi — > on the left-hand side of this equation; it includes all 
pole terms plus a non-singular contribution that restores gauge invariance. Only these Born terms contribute to the 
low-energy limit of Compton scattering because they are of zeroth order in the photon momenta; we therefore find 

Hm T»» = Hm (T^ LE + T& po J = Hm Wx (82) 

and thus end up with the correct Born-Thomson limit for the Compton scattering amplitude if we use the full 
two-photon irreducible vertex Y^ v : 

There is also an alternative way to demonstrate the correct low-energy behaviour of the full two-photon irreducible 
vertex r MI/ . We recall that we have only used results from differentiations of amputated (point-form) operators such 
as J' 1 and Y^. The pole terms can also be found if we define in analogy to J M as a twofold differentiation of the 
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free boson propagator Ap(P) with respect to the momentum P followed by an amputation of the resulting vertex in 
the limit q — > (see e.g. P0| ): 

f" v ■= limT"" = A F * \{-Qd v ) (— Q9 M ) A F ] A f 1 = J^AfJ 1 ' + J^AfJ" - Q 2 ■ 2 ff ^ (83) 

Inserting the explicit expressions for J M and A p , we find the low-energy limit of the Born term in eq. ( |75| ) including 
all pole contributions. This procedure is equivalent to a differentiation of the full expression — xT^X m ec l- (|7^)i in 
this way, we do not only take into account the derivative of the amputated (point-form) vertex but also include 
terms proportional to d v \ and d v \ which will then re-introduce the pole terms that are also present in T^ v . With 
the differential identities for d v T^ and for the Bethe-Salpeter amplitudes in eqs. @, @ and (§|), and fr om the 
vertex' re-definition in eq. ([72]), we recover 

- (-OP) x^x = Q({d v x)^x + x{d v T^) x + x^(d v x)) 

= iX^ v X ■ 

For the relevance of the formal differentiation of the Bethe-Salpeter amplitudes X and x which are properly defined 



only on mass shell, see the remarks in section IV E. It is crucial for obtaining this correct result that the Ward identity 
of the order 0(e 2 ) in eq. ( |5^ ) includes subtraction terms that just cancel the contributions from d v x and d v %. In 
other words: the pole terms are re-introduced by the contributions coming from the Bethe-Salpeter amplitudes in 
an external field, see fig. (|l^). In this sense we can state that a consistent differentiation prescription with regard to 
point-form parts (i.e. J M and T M ) and non-amputated parts (like Bethe-Salpeter amplitudes and boson propagators 
that are cut after taking the derivatives) leads to the correct Born-Thomson limit. Both ways yield, as shown above, 

S% T " V = ft, ( T Po- + T NoPole) = ixT^X ■ (84) 

However, we shall note here that we prefer the first derivation which can be formulated without using terms like d u x 
and d v x which are not strictly defined for on-shell quantities like Bethe-Salpeter amplitudes; note that the relations 
in eqs. ( |6l| ) and (^3|) have been found on quite formal grounds. 

The Born-Thomson limit can also be derived if we only include re-scattering terms due to T matrix contributions but 
neglect one-photon or two-photon irreducible interaction kernels. This can be seen by applying the same procedure 
as before to the lowest order current matrix element in eq. ( f74| ) : 

^NoPole = -2 2 ' = -OPJ" = ~X( - Qd v T») X (85) 

= - r%G r» - r^G r(;) x 

here, we have only considered the differentiation of the amputated (point-form) operators and used that G = G Q — iT . 
We thus find — > ^F^" as the modified tensor defined in eqs. (pi]), ( p9| ) and ( [72] ) as well as pole terms that can 
be re-arranged to the left-hand side. For — » 0, the pole terms in the Green's function G dominate all other 
contributions; analogously to the derivation of the low-energy limit of the full two-photon vertex, we are led to 



H (* r o X ■ ■ XV U X + xT'oX ■ ■ xtfZx + r^ P0LE ) = Hm ix^X ■ (86) 



Since the one-photon irreducible interaction current (which is neglected in this approximation) does not contribute 
in the static limit for interaction kernels K that are independent of the total momentum, see also eq. (f74|), we find 
the correct low-energy behaviour for the two-photon vertex T^: 

Hm T<- = fim (T p ^ LE + T^ Po J = hm i x t^ X ■ (87) 

This leads to the Born-Thomson limit for the Compton amplitude although we neglected contributions from n— photon 
irreducible interaction kernels in T'"' := ixT^x'- 



2 



lim 3V = - JL Hm e^e 2v = (ft • ft) 
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Let us finally note that we cannot expect to recover the correct low-energy limit if we only consider the lowest order 
contribution to Compton scattering in Tq 1 := ixPoX' 



1 Q 2 



This is obvious from the considerations given above: if we do not include pole contributions then we will not be able 
to produce the full gauge invariant Born term, see eq. ([75]). However, in contrast to the — Q ■ 2g^ v part, these pole 
terms do not contribute to the Compton amplitude T; but for the correct reproduction of the limit — %- (ei ■ e%) we 
must at least include the full non-singular tensor I^oPolb anc ^ n0 * on ^ ^ s l° wes t-order part, see eq. (|80|). 
As an example, we will finally write down explicitly this lowest-order Compton scattering tensor Tq" — ix^o X f° r 
a qq system with finite photon momenta q±, q2 0: 



\qq 



ie\e 2 



ie\e 2 



ie\ 



ie\ 



ie% 



d 4 p 



(2tt)< 

d 4 p 

d 4 p 



X(P,P + m<li +»7i92) 

4 x(P,p-vm -mq2) 



(2tt) 4 
d 4 p 

(27T) 4 

(2ir) 4 
d 4 p 



x(P,p + mn - V2q2) 

x(P,p-vm +vi<i2) 
4 x(P,p-vm +V1Q2) 



r®Y)x(p,p) 
Y®r)x{p,p) 

'fSpfaP + P + 9i)7 M ® Sp^-mP + pj) X(P,P) 
rS F (r ]1 P + p-q 2 )Y^S F 1 (-V2P + p))x(Pp) 

S F \ m P + p)(g> YS F {~mP +P~ qih") x(p,p) 

s F \mP + p) ® rs F {- m p + P + q2 )Y) x(p, p) 



Again, the correct charge factors are found for the different terms contributing to Tq 1 "; note that the anti-quark with 
momentum — p 2 = —P/2 +p has the charge — e 2 . Six diagrams emerge from this lowest order two-photon vertex; if 
T matrix contributions are taken into account and interaction currents are still neglected, eight additional terms will 
have to be considered for a qq system. 



VII. SUMMARY AND CONCLUSIONS 



We have derived the second order Green's function G 77 for Compton scattering off a bound state in a scheme that the 
authors of ref. ji]-[§ (where it has been introduced for processes in first order of the electromagnetic coupling) have 
labelled "gauging a hadronic system". The resulting full two-photon vertex T 77 = G _1 G 77 G _1 includes the lowest 
order contribution Tq 7 and an explicit two-photon irreducible interaction kernel .ft" 77 but also re-scattering terms via 
an intermediate full Green's function r 7 Gr 7 (with G = Go — iT) as well as subtraction terms r^Gor^. Note that 
the singularities in r 7 Gr 7 due to the poles in the Green's function for vanishing photon momenta are responsible for 
the correct low-energy limit of the Compton scattering tensor including the pole contributions present in the Born 
terms. 

By imposing that the amplitude of order 0(e 2 ) obeys the related Ward-Takahasi identity, we found a constraint for 
the two-photon irreducible interaction kernel similar to the one in first order (see also JlM). We also dicussed the 
Ward-Takahashi identities for the "gauged" Bethe-Salpeter amplitudes x M and x^ h ' ■ Differential identities for Green's 
functions, the vertices and and the n-photon irreducible interaction kernels were derived. By using rather 
formal arguments, we found similar identities for the Bethe-Salpeter amplitudes; however, we stressed that these are 
strictly defined only on mass shell, i.e. for P 2 — M 2 . 

As we checked the gauge invariance condition for a matrix clement proportional to xr 77 x, we found that it is (as 
it should be) automatically satisfied if the two-photon vertex obeys the correct Ward-Takahashi identity. For an 
approximation in which wc neglected all n-photon irreducible interaction kernels, we found that gauge invariance for 
the resulting two-photon vertex T 77 is preserved only if the Bethe-Salpeter kernel K does not depend on the bound 
state's four-momentum P and is of local type, i.e. K(P; {h}, {fc^}) = K({ki — k^}) where {hi} denote the relative 
momenta of the constituents. This is the same condition as for the related approximation in first order if we neglect 
the one-photon irreducible kernel, i.e. r 7 = T 7 = Tg. 
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A similar result could be derived with regard to the low-energy limits: they are properly described by the full vertices 
r 7 and r 77 . If we only neglect contributions from K 1 and K 11 , then we also meet the correct limits for vanishing 
photon momenta. However, we do not find the correct Born-Thomson limit for the Compton scattering amplitude 
if we only include the lowest order irreducible vertex T 7 , 7 . Therefore we conclude that the corresponding vertex to 
the lowest order vertex T 7 = r 7 , at 0(e) is in second order of the electromagnetic coupling not Tq 7 but merely 
T 77 = Tq 7 — iI^Tr 7 , as a vertex that also includes re-scattering terms due to T matrix contributions. 
Let us finally note that the scheme presented in this article can also be applied to processes like 77 — > AA, i.e. the 
production of bound state pairs {e.g. it + tt~ or 7r°7r°) in photon-photon collisions. The corresponding identities could 
as well be obtained by introducing new momentum variables (72 - * ~Q2 and P — > —P that satisfy the four-momentum 
conservation q\ + q% = P + P'. Furthermore, we want to stress that our results could also be of interest for Compton 
scattering off an electromagnetic bound state such as positronium since we didn't specify the interaction kernel in the 
Bethe-Salpeter equation, and in this respect our considerations are quite general. 
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APPENDIX A: COORDINATES FOR A COMPOSITE SYSTEM 



Let us consider a system of n fermions and n anti-fermions (N = n + n) where the position of each constituent is 
described by the coordinate Xi (z = 1,2,..., N). Now we introduce new coordinates by 



X := rjixi 



VnXn 



(Al) 



ri :— xi - x 2 

r-z := ~(x\ + x 2 j - x 3 



1 



rjv-i 



-(xi + X 2 + + ■ ■ ■ + XN-l) - Xn 



N —1' 

The coefficients rji have no direct geometrical meaning; the only constraint is that they have to sum up to unity: 

rji + T)2 + ■ ■ ■ + r) N = 1 . 

Choosing them as T]i — rrii/M with M = X]i m i gi yes X the meaning of a center-of-mass coordinate; however, a 
more common choice is simply rji = 1/N. Another possibility is to interpret X as a center-of-charge coordinate by 
adopting r]i = e.i/0. where ej is the physical charge of the i-th constituent and Q = X)i e * i s the total charge; this 
choice will be applied in the last section of this article. Note that the fixing rji = ei/Q is formally also possible for 
neutral bound states since even for Q — the rela tion ^ . rji = 1 is satisfied. 
In a compact form, the transformation in eq. (Al) can be written as a matrix equation like 



with a. 



( X \ 



( Vl V2 

1 -1 



and A x 



(3 X 



V3 





V N-l N-l N-l 



X2 

\x N J 

VN \ 






- 1 / 



(A2) 



Now consider the canonical conjugated momenta pi,P2, ■ ■ ■ ,Pn and P,ki, . . . , fcjv-i, respectively, where P is the total 



momentum P = p\ 
well: 



P2 + ■ ■ ■ + Pn- We can give a matrix equation for this transformation in momentum space as 



(A3) 



Obviously, the z-th element in the first column of 1 gives the dependency of the momentum pi of the total 
momentum P. 

The transformations given in eqs. ( |A2| ) and ( [A3] ) satisfy the following relations concerning their Jacobi determinant 
and their scalar product: 





( P \ 




(n\ 








P2 


(*p = Ap/3 p with a p = 














{p N ) 



d(X,n, . . . ,r N -%) 



d(xi,x 2 , ■ • ■ ,Xn) 



d(P,h 



<9(pi,p 2 , ■ ■ ■ ,Pn) 



and 



PX + r\kx + . . . r N _ik N _i = p 1 x 1 + p 2 x 2 + ■ ■ ■ + PnXn ■ (A4) 
The last equation gives the transformation matrix A p in momentum space in terms of the matrix A x in position space: 
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a\a x = (3 l p l3 x a\k x $ x = (A p l a p ) t p x = a* (A; 1 )'^ . 

This yields A" 1 = A^, so that we find the dependency of the momentum of the i— th constituent of the total momentum 
to be 

Pi =T)iP + f(k u ...,k N -i) (A5) 

where f(ki, ■ ■ ■ , fcjv— l) is a linear function that can be derived from A" 1 = A^j,. The constituents of a qq meson for 
instance have the momenta 

Pi=7/iP + p and —p 2 = —iq%P + p (A6) 
with total momentum P = pi + P2 and relative momentum p := k\ = ?/2Pi — t?iP2- 

APPENDIX B: BETHE-SALPETER AMPLITUDES IN SECOND ORDER 

In this appendix, we want to derive the Bethe-Salpeter amplitude in second order of an external field. We start with 
eq. ( fjl| ) and write down the expansion in eq. (fjj) for = Go, K, x- 

X - iex 1 - e V 7 + ■ ■ • = -i(G - ieG 7 - e 2 G 77 + . . .) [K - ieK~< - e 2 K^ + . . .) ( X ~ ie X ~< - e V 7 + . . .) . 

In the order C(e 2 ) of the electromagnetic coupling, we find with the first order relation x 7 = Gr 7 x 

X 77 = - t G^Kx - tGoK^x ~ iG^Kx 11 - iG^x ~ iG&Kx 1 ~ zG ^ 7 x 7 

= (l + iG K)-\ (-iGoT^GaKx-iGoK^x-iGoTlGoK^x-iGoTl G KG T^ x - iG K^GT~<x 

=GG- 1 = l X =»(G-G ) 

= g(t 77 - iK^ - r 7 G (r 7 + iK~<) + (r 7 - ix 7 )Gr 7 ) x 
= G(r 77 - iK^ + r 7 Gr 7 - r 7 G r 7 )x 



where we have used the fundamental equation for the Green's function, the Bethe-Salpeter equation and some relations 



introduced in section [II. With r 77 := T 77 - iiH 7 + r 7 Gr 7 - r 7 G r 7 , we finally find the expression in eq. (§3) for 



the second order Bethe-Salpeter amplitude x 77 = Gr 77 x. 
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